Exponentially fast Thinning of Nanoscale Films by Turbulent Mixing 
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Films are nanoscopic elements of foams, emulsions and suspensions, and form a paradigm for 
nanochannel transport that eventually tests the limits of hydrodynamic descriptions. Here, we study 
the collapse of a freestanding film to its equilibrium. The generation of nanoscale films usually is a 
slow linear process; using thermal forcing we find unprecedented dynamics with exponentially fast 
thinning. The complex interplay of thermal convection, interface and gravitational forces yields 
optimal turbulent mixing and transport. Domains of collapsed film are generated, elongated and 
convected in a beautiful display of chaotic mixing. With a timescale analysis we identify mixing as 
the dominant dynamical process responsible for exponential thinning. 



PACS numbers: 47.61.-k,47.57.Bc,82.70.Rr, 47.51. +a 

In nature, thin liquid films form the interfaces in foams, 
emulsions, or they sheath solid particles in a liquid en- 
vironment. Currently, interest in thin films is driven by 
science, but also by technological demands, since they 
are the basic structural elements of dispersed systems, 
separating gas bubbles in foams, droplets in emulsions 
and particles in suspensions. 

Thin film dynamics is governed by gravitational^ cap- 
illary or interfacial forces, inducing the existence of 
the disjoining pressure. The latter combines long- and 
short-range molecular forces: electrostatic, van der Waals 
(vdW), and steric forces [1, 2 and strongly depends on 
the distance between the interacting surfaces. Whereas 
films on substrates are established in industry and re- 
search [3 , freestanding thin liquid films still provide a 
challenge in experiments and theory alike. Consequently, 
the study of foam films is of central interest e.g. [4]-[7] 
We contribute by presenting a novel approach to thin- 
ning of vertically oriented, freestanding, thermally forced, 
nonequilibrium foam films. 

As a result of the force balance two stable equilibria 
may occur, depending on the chemical composition of 
bulk solution and chosen surface active agents (surfac- 
tants): Common Black Films with a thickness of more 
than 10 nm are formed when electrostatic interactions 
balance the dominant van der Waals foce, and, of course, 
gravity and capillarity [U |8]; Newton Black Films are 
stable with a thickness of less than 10 nm, due to repul- 
sive short range steric forces O |10] . Vertically oriented 
foam films are initially thick (500 - 5000 nm), but in 
the thinner regionsbut after some time, an area of black 
film (BF), which can be the Common or Newton type 
appears uppermost, with a sharp horizontal boundary 
towards the thick film, see Fig. [l] 'E'. As the foam film 
drains, this boundary moves downward. Since the menis- 
cus at the boundaries is typically much thicker, the main 
transport happens there [11 . The final BF is thermody- 



namically stable. 

The downward motion of the BF front is at constant 
velocity, limited by the thin film bulk velocity which fol- 
lows a Poiseuille flow [12]. The drainage, steady-state 
profile and thermal effects have been predicted earlier , 
experiments with evaporation and comparatively uncon- 
trolled chemistry [13, 14 suggest an exponential decrease 
in thinning speed. 

Here, a foam film is brought far into nonequilibrium 
by pointlike thermal forcing, cf. Fig. [l} b. The addi- 
tional force counteracts gravitation and capillarity, and 
changes the dynamics from linear to exponential. After 
a transient period, a global convective flow with com- 
plex dynamic behavior is established, marked by two- 
dimensional turbulent mixing combined with surface in- 
stabilities. Image recording allows for a characterization 
of the flow field, global mixing properties, and eventually 
the thinning law, as the film progresses towards a global 
BF state. As it turns out, the observed behavior is only 
possible by a complex interplay of various dynamical pro- 
cesses within the nanoscale film. 

Experiment The experimental setup consists of a ver- 
tical rectangular aluminium frame with rounded corners, 
45 X 20 mm, enclosed by an atmosphere-preserving cell 
with a glass window for video recording, cf. Fig. [l] more 
details and video, see [15]. Both Common and Newton 
Black Film are optically transparent, and show only a 
subtle difference in reflectivity which is not resolved by 
our setup. For the set conditions it is expected that a 
Newton Black Film forms, however it will be referred 
to with the general term black film (BF). Thermal forc- 
ing is effected by inserting a cooled copper needle at the 
film center (T = — 169°C). The needle enters the cell 
through a fitting hole, centred on the thin film area. Am- 
bient temperature was constant at 20 °C, with a Rayleigh 
number Ra ~ 10^, typical for turbulence. 

The solution from which the liquid film was drawn con- 
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FIG. 1: The (left and right) snapshots show the reflection 
of a diffuse light source off the foam film surface. The color, 
as seen by the naked eye, repeat from top to bottom with the 
increasing film thickness, a Foam film in aluminium frame - 
thinning without thermal forcing, b Perspective view . A: 
Light source, B: diffusor screen, C: video camera, D: casing 
to prevent evaporation, with plexiglass window, E: BF area 
F: Cooling copper needle at — 169°C, G: Soap reservoir, if; 
Aluminium frame The convective pattern, established with 
thermal forcing is sketched, along with the BF and the thick, 
wedge-like region. The latter decreases in the course of time, 
c Foam film in aluminium frame - thinning with thermal forc- 
ing. The convection rolls are in the lower half of the frame, 
above the cold needle a stably stratified region has formed. 



sists of the surfactant n-dodecyl-/3-maltoside {^^—Cx^G^)^ 
prepared with filtered deionized water and stabilized with 
25 %voi glycerin [16 . The entire film area is observed op- 
tically by a conventional video camera [15 . 

The interference of incident and reflected light yields a 
striped interference pattern, which can be used to infer 
the film thickness. Each color cycle corresponds to the 
multiples n of the smallest negative interference condition 
n\ = 4/i. In principle, reflection depends on temperature 
via the refraction index, here it is assumed constant. The 
velocity is measured by CIV, color imaging velocimetry 
[17]. The main concern is the thinning speed, so extra 
care is placed on its exact determination: the color tables 
and contrast of the video images were tuned to achieve 
pixel-level accuracy of the position of the interface be- 
tween BF and thick film. 

Results Without thermal forcing, the initially thick 
film (micron scale) is drained by gravitation and capil- 
larity, thereby evolving towards its equilibrium BF thick- 
ness. This thinning is limited by Poiseuille flow [131 HH] 
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with the velocity given hy v In the bulk, 

h - 10-^ - 10-^ m, and v - 10"^ - 10"^^ ^/s. On the 
vertical boundaries, the so-called Plateau border with a 
meniscus of h ^ 10~^ — 10~^ m allows for a faster trans- 
port at 1' 10~^ — 10~^ ^/s. The film in this frame 
reaches equilibrium, i.e. complete BF phase after approx- 
imately 5 h. A corresponding series of snapshots is given 
in Figj2j a; the dominant transport can be observed by a 
Kelvin-Helmholtz-like instability at the Plateau border. 
We clearly observe a linear behavior with constant ve- 



locity Vreguiar — 18 /^"^/s, in accordaucc with [12]. The 
very first quantitative thinning experiments [M] reported 
an exponential slowing down. This can most likely be at- 
tributed to uncontrolled evaporation from the film since 
the atmosphere was not saturated, already anticipated 
by Mvsels[T4]. For the aspect ratio here (5:3), the main 
flow is through the Plateau border, as explained above. 

With thermal forcing ((A) in Fig. [2] and Fig. [3|, the 
film is cooled down rapidly in the center and a small, 
frozen region appears surrounding the needle. Thermal 
energy is exchanged at its boundaries and heavy convec- 
tion starts. A weakly turbulent flow is rapidly established 
with 2 symmetric, dominant eddies, cf. Fig.[l] c. As a re- 
sult of this convection, two different regimes are observed: 
a first regime with faster, but still linear thinning, and a 
second regime with exponential thinning. 

In the linear regime (150 s-277s) the velocity of the BF 
front is still constant, however the magnitude is increased 
by a factor of 3 to vi =66^^/^ (cf. Fig. [3) blue line). The 
enhancement is explained as follows: 

Fluid is transported by convection vertically from top 
to bottom, however the outflow into the reservoir is 
slower than convection. Necessarily a vertical upflow at 
the left and right boundary is generated. This in turn 
causes a suction of fluid from the bottom. This results in 
faster uniform motion of the BF front, whose speed is now 
determined by the net deposition rate of fluid through 
the bottom meniscus and subsequent outflow from the 
same, and is no longer limited by the flow through the 
vertical menisci. The latter is enhanced by the more con- 
stant thickness of the film (with height) resulting in the 
increased velocity. 

The second, exponential regime starts at t 277 s (B 
in Fig. [2|, its dynamics is explained as follows: The BF 
spots arise spontaneously due to thermal and mechanical 
fluctuations, and are subsequently transported around by 
the convection. A BF spot is basically a domain of BF 
(5-20 nm) surrounded by much thicker film (100-1000 
nm). This BF spot can only grow by transporting liquid 
away from the spot through the surrounding film. Hence, 
this must take place under Poiseuille flow conditions, due 
to the nanoscale thickness of the film. Du to gravity, a 
BF spot lifts up, like a bubble. 

Due to convection, the growth of BF spots is affected 
by the shear stress at their interface. Beyond a certain 
threshold of size and speed (approx. 285 s, see also [15] ) 
the BF spots are stretched, leaving behind a thin BF fil- 
ament. These filaments are folded and seed i) new spots 
acting as multipliers, ii) "droplets" of thick film within 
the BF spots. Consequently two additional flows are es- 
tablished: BF spots are lifted to the top to join the BF 
area, droplets descend due to gravity and eventually flow 
out into the bottom reservoir. The droplets appear in 
a grey color since their reflection is diffuse, (cf. Fig. [2j 
285 s). These processes can only take place when the BF 
spots (droplets) are light (heavy) enough to escape the 
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FIG. 2: Image sequences of foam film thinning. Thinning behaviour of a foam film a without thermal forcing, and b with 
thermal forcing. While the undisturbed thinning (top) evolves slowly, the thermally driven flow (bottom) exhibits convection, 
with rapid transport and mixing. A: start of convection, B: onset of exponential thinning. One clearly recognizes the black 
film spots which are stretched and folded, such that a rapid conversion the film evolves to the black film equilibrium phase. 




FIG. 3: Temporal evolution of black film area: linear 
versus exponential thinning, grey triangles: without 
thermal forcing, the thinning is perfectly linear in time, blue 
circles: with thermal forcing, one observes three different 
states: for times t < 150s (transition at A), no thermal force 
is applied and the convection has not yet started. For times 
150s < t < 280s, the convection is established, but BF spots 
are not stretched into filaments, consequently thinning is lin- 
ear but faster. Finally, at t = 280 s the transition to exponen- 
tially fast thinning can be observed, marked by B. BF spots 
and turbulent mixing are acting inside the convection zone. 
Solid lines are linear fits. The inset demonstrates the expo- 
nential behavior by a high-quality fit for 1 — A = exp(— t/r), 
with a timescale r =4.76 s. 



convection roll. The transport of the droplets is much 
faster than Poiseuille flow permits, and the exact motion 
is therefore difficult to describe analytically. 

Let us comparing the dominant timescales for thinning 
with and without thermal forcing. Essentially, the thick- 
ness varies due to i) density variation according to the lo- 



cal temperature; this effect is of the order 
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ii) advection of h with a ratio 5/500 (final vs. initial thick- 
ness). 

This thinning process is multiplicative and effectively 



obeys an exponential decrease of the transient, thick 
phase: A = 7(^0 — A), with A the area covered by 
thick film, a time constant r = I/7 and a thinning speed 
V A/Ly with Ly the lateral extension of the film. Space 
directions are as indicated in Fig. [l] b. The thinning 
timescale is then r which is measured to be r =4.76 s, a 
thinning process which is 500 times faster than usual! 

In general, the timescale depends on geometry and flow 
properties and enters parametrically into r{v, L^^ Ly). 
To obtain an estimate, consider V = J^^dxdy h{x,y), 
where Aq = - Ly is the total area. The area-averaged 
growth speed is given in the most unstable, varicose mode 
12. 
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dxdy h = — / dxdy h{x^y)S/ 2u{x^y) (1) 
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with V2 the two dimensional divergence and u the corre- 
sponding velocity field. Note that i) h is the convective 
time derivative, i.e. the local growth and ii) the flow is 
highly compressible in 2D, justifying the above equation. 
On average, as denoted by brackets, the local thickness 
becomes h = —h {V2U) ,, which is coupled to a compli- 
cated equation for the velocity. 

The flow is dominated by two local structures, the BF 
spots and the BF channels. Both expand and thus {V2U) 
is larger than zero, i.e. the thickness h decreases exponen- 
tially. In the absence of a restoring force, this would lead 
to rupture; in this case the divergence of the velocity con- 
verges to zero for the thickness going to its equilibrium 
BF value. However, if the thickness decreases in one area, 
there must be a corresponding increase in nearby regions, 
leading to "rims", similar to the ones found in evapora- 
tive dewetting [19 . They are also subject to gravitation, 
and fluid flows to the bottom. The time scale is now set 
by the convection in contrast to the Poiseuille flow for 
the unforced case. 
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This can be tested by calculating the time scales of the 
involved processes: thermal forcing, gravitation, steric 
and electrical forces (including Van der Waals), cap- 
illarity, mechanical strain (dissipation). To estimate 
the important effects, we calculate Rayleigh-, Capillary-, 
Marangoni-, Biot-, and Prandtl numbers (i?a, Ca, Ma, 
5, Pr), and the dimensionless Hamaker constant A. 
This leads to the following values: Ra ~ 10^, Pr ~ 7, 
Ma - 10^, P - 10-^, A - 10-"^, from which we de- 
duce that the dominant force is thermal. The convec- 
tive time scale is calculated by the characteristic turnover 
time 5 s. Eventually, the time scale of the 

strain, Tstrain ^ 0.4 5, responsible for the local genera- 
tion of BF spots. It is smaller than the convective one, 
consequently convection rules the exponential decrease of 
thick film. The measured time constant r = 4.76 coin- 
cides perfectly with the convective time scale, given the 
roughness of our assumptions and the complexity of the 
process. 

To produce a local BF spot, the surfactants must fluc- 
tuate very rapidly such that a hole cannot be filled im- 
mediately from the bulk excess concentration. Stretching 
and folding: from the local strain rate Tstrain^ the mix- 
ing efficiency, [20], E = 0.075, is calculated. Folding is 
related to the intersection of stable and unstable mani- 
folds, which in our experiment are time dependent. The 
time-averaged flow has two hyperbolic points, which are 
heteroclinically connected, providing a basis for a horse- 
shoe providing the folding. 

Discussion We have demonstrated the dramatic ac- 
celeration of the thinning of a foam film by thermal forc- 
ing. We use novel experimental techniques to control the 
properties of the thin liquid film in a precise manner, i.e. 
a measurement chamber that inhibits evaporation from 
the film surface. The analysis of the fluid flow is local and 
global, by tracking BF spots to measure local stretching, 
folding and global thinning rate. 

The thermal driving leads to convection, which changes 
the dynamics of thinning from linear to exponential. This 
is not only substantially faster, but qualitatively differ- 
ent from the classical, thermally homogeneous, thinning 
process. This result can be explained: the growth of 
black film area is driven by the pressure gradient at the 
interface with the non-equilibrium phase. Thus the in- 
terface length determines the overall thinning rate. By 
advection and stretching of black film spots, trails are 
left behind which increase multiplicatively the interface 
length resulting in the exponential thinning behaviour. 
Efficient transport and filamentation are guaranteed by 
the mixing properties of the flow. 

We plan to widen our work to a Rayleigh-Benard setup 
where the influence of interfacial forces can be studied, 
or, more technically, in a controlled mixing of substances 
in a quasi 2D setup. Further, with chemical control one 
can tune the film thickness and thereby observe the de- 



pendence of thermodynamical and hydrodynamical prop- 
erties on the thickness, eventually testing the limits of 
continuum description. A comparison with quantum cal- 
culations will be a formidable task for future research. 
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